Bulletin of Electrical Engineering and Informatics 
Vol. 9, No. 4, August 2020, pp. 1685~1693 
ISSN: 2302-9285, DOI: 10.1159 1/eei.v9i4. 1886 O 1685 


Optimization of distribution network configuration with multi 
objective function based on improved cuckoo search algorithm 


Thuan Thanh Nguyen 
Faculty of Electrical Engineering Technology, Industrial University of Ho Chi Minh City, Vietnam 











Article Info ABSTRACT 

Article history: This paper proposes an improved cuckoo search (ICSA) for solving 
. the distribution network reconfiguration (NR) problem with multi-objective 

Received Nov 2, 2019 function. The membership functions are considered consisting of minimizing 

Revised Feb 6, 2020 of power loss, load balancing among branches and among the feeders, node 

Accepted Apr 1, 2020 voltage deviation and switching operation numbers. ICSA is developed 


from the original CSA with adding the local search mechanism for exploiting 





around the current best solution. The effectiveness of the ICSA has validated 
Keywords: on the 70-node and the 83-node practical systems. The obtained results have 
been compared to those from runner root algorithm (RRA) and other 
methods in the literature. The obtained results demonstrate that ICSA 
Improved cuckoo search has high ability for searching the optimal solution with higher successful rate 
Multi-objective and better quality of obtained solution as well as smaller iterations compared 
Network reconfiguration to RRA and other methods. Therefore, ICSA is a reliable method for 
Power loss reduction the multi-objective NR problems. 


Electric distribution network 


This is an open access article under the CC BY-SA license. 





Corresponding Author: 


Thuan Thanh Nguyen, 

Faculty of Electrical Engineering Technology, 

Industrial University of Ho Chi Minh City, 

No. 12 Nguyen Van Bao, Ward 4, Go Vap District, Ho Chi Minh City, Vietnam. 
Email: nguyenthanhthuan @iuh.edu.vn 








1. INTRODUCTION 

Transferring a part of loads from the heavy branches to the light branches is not only reduction 
of power loss but also gains a lot of advantages such as improvement of voltage profile, balance of load 
among branches and the feeders and enhancing the reliability of the electric distribution systems (EDS), etc. 
Transferring loads among branches in the EDS is implemented by changing the status of switches located in 
each branch. This processes is called the network reconfiguration (NR) problem. The NR _ problem 
is a nonlinear problem with constraints and a lot of local extremes. Therefore, studying the methods 
for solving the NR problem have attracted many concerns of the researches in more than for decades. 
From the first proposed by Merlin and Back [1], the NR problem has been solved by a lot of methods with 
different objectives. 

In [2-4], the methods based on the heuristic search techniques have proposed for the NR problem 
with multi objective functions such as reduction of power losses, enhancement of the node voltage 
and increase of reliability index [2], reduction of power loss and the maximum current of the EDS [3] 
and optimizing the bus voltage deviations, the branch current loading, the phase current deviation and power 
loss of the EDS [4]. In comparisons with the heuristic methods, the number of methods based on 
the meta-heuristic techniques is much greater than due to a lot of advantages such as easy to implement 
and apply for different kind of the objective functions and the better obtained results. Therefore, more 
and more the meta-heuristic method are proposed for the multi objective NR problem. In [5], discrete 
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particles swarm optimization (PSO) algorithm has proposed for finding the optimal configuration of the 33 
and 69 nodes systems to satisfy the objective functions consisting of minimization of power loss, the number 
of switching operations and deviations of bus voltages. In [6], PSO is applied for the multi-objective 
functions consisting of reliability indexes and reduction of power loss. In [7], Binary PSO is used for finding 
the optimal configuration to decrease power loss and the voltage stability index. Harmony search algorithm 
is proposed for the NR problem for minimizing of real power losses and improvement of node voltage. In [8], 
adaptive genetic algorithm is used for determining the best configuration of the 69 and 136 nodes systems 
with the objective functions of reduction of power loss, node voltage deviation, branch current loading index, 
the switching numbers. In [9], the multi objective function consisting of minimizing the power losses and the 
reliability indices is solved by the genetic algorithm (GA). In [10], GA is also is used to solve the NR 
problem for reducing power loss and costs of maintenance and operation of switches. Micro GA [11] is also 
successful applied for finding the optimal configuration to minimize power losses and enhance a reliability 
index of the EDS. In [12], honey bee mating optimization is presented to reduce power loss, the switching 
numbers and the voltage deviation for the 33 and 86 bus systems. In [13], The objective functions such as 
system average interruption frequency index, average energy not supplied and active power losses are 
optimized by bat algorithm. In [14], grey wolf optimizer algorithm is proposed for the combined problem of 
network reconfiguration and STATCOM placement to reduce power loss and load balance index. The 
shuffled frog-leaping [15] is used for the problem of network reconfiguration with reactive power 
compensators. In [16], evolutionary programming, firefly, PSO and gravitational search have been successful 
applied for the multi objective reconfiguration to reduce power loss, enhance voltage profile and maximize 
the distribution generator output. In [17], the network reconfiguration considering reactive power dispatches 
of distributed generations (DG) is solved by non-dominated sorting PSO to reduce energy wastage from DG, 
power loss and voltage deviation. Reduction power loss, maximum voltage deviation, load balancing index 
and switching numbers of the 33 and 84 nodes system are considered by invasive weed optimization [18]. 
Reduction of the power loss and maintain of a continuous power supply to the customers are also considered 
as the multi objective function of the NR problem an solved by modified ant colony algorithms [19]. From 
above studies, it can be seen that with the ability for providing a good solution to NR problem, more and 
more recent meta-heuristic algorithms will be applied for the NR problem. 

The cuckoo search algorithm (CSA) is a recent meta-heuristic inspired from the cuckoo birds’ 
reproduction [20]. From the first proposed in 2009, CSA have been applied for a lot of problems related 
to the power system such as the hydrothermal scheduling problem [21], the problem of heat and 
power economic dispatch [22], the problem of placement of distributed generators and capacitors in 
the EDS [23, 24]. For solving the NR problem, CSA has been successful applied for the problem to reduce 
power loss and improve voltage profile [25]. In this study, CSA outperforms PSO. In [26], CSA has been 
applied for the NR problem considering placement of DG. In addition, to enhance the efficient of CSA for 
the NR problem, the improved CSA (ICSA) has been proposed for solving the NR problem to minimize 
the power loss [27]. In this study, ICSA outperforms to CSA as well as many improved versions of 
CSA. However, the efficient of ICSA for the NR problem with other objective functions is still 
a considerable question. 

This paper presents an NR method for multi-objective based on ICSA. The objective functions 
considered are reduction of power loss, enhancing of load balancing among the branches and among 
the feeders, node voltage deviation and the number of switching operations. The effectiveness of the ICSA 
for the multi-objective NR problem is evaluated on the 70-node and the practical 83-node systems. 
The highlighting features of the work can be summarized as follows: (1) The ICSA is proposed for 
the multi-objective NR problem with five above membership objective functions; (2) The selection of 
the final solution based on max-min method is used to compromise the membership objective functions; (3) 
The performance of the ICSA is effectively evaluated on the complex EDS consisting of the 70-node and the 
practical 83-node systems. 


2. THE PROBLEM FORMULATION OF THE MULTI OBJECTIVE NR 

There are five membership functions are considered in the paper consisting of minimizing power 
loss, improving load balancing among the branches, load balancing among the feeders, improving node 
voltages and reducing the number of switching operations. Each membership function is determined as 
follows: 


2.1. Reduction of power loss (AP) 
Power loss reduction is the most part of the NR problem which is determined as follows: 
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Where Np, is the branch number of the EDS. Ri, Vi, P; and Q; are the resistance, the ending voltage, the active 
and reactive power flow of the ith branch, respectively. 


2.2. Balance of loads among the branches (LBI) 

Balance of loads among the branches helps the load carrying level of each branch in the EDS 
is approximately equal. Due to the different rate current of each branch, to balance the load carrying level of 
each branch, LBI is determined as follows [18, 28]: 


I I Ij In 
LBI = var es 
I1,max 12,max Timax Inmax 











(2) 


Where var called variance, J; and J,max are the current and the rate current of the branch ith, respectively. 
2.3. Balance of loads among the feeders (LBF) 


LBF is considered to support the EDS transfer loads from the heavy feeders to light feeders, 
it is determined as follows: 


EBF = vor[lpideoidey es leyel (3) 
Where J; ; is a current of the ith feeder, NF is the total number of feeders. 
2.4. Node voltage deviation 

To enhance power quality of the system, the node voltage deviation (AV) is also considered to an 
objective of the NR problem. AV is calculated as follows: 

AV = Veer — Vinin (4) 
Where V,-¢¢ and Vinin are the feeders voltage and the minimum node voltage in per unit of the EDS. 
2.5. Number of switching operations (NSW) 

For decreasing the operating cost, NSW is also considered as one of the objective function. 
Which is determined as follows [18]: 

NSW = YT|Soi — Si| (5) 
Where, So; and S; are status of the switch ith before and after performing reconfiguration. 
2.6. The objective function 


The final compromise solution (FCS) for five membership functions (MF) is determined by 
the max-min method [18, 29, 30] as follows: 


1, Rios re" 
MF,(X) = ae Fmin < F(X) < Fm (6) 
0, BO) 2 


In which, the maximum value of the membership functions Fz"** are calculated from the initial 


configuration because the purposed of the NR problem is finding the configuration having the better 
objective function value than the initial configuration. The F7”" of the membership functions are determined 
by solving the NR problem with single objective function. 

The degree of overall satisfaction for each candidate solution is the minimum of all of 
the membership functions. Then, due to the membership functions value of less than unity, the fuzzy decision 


for the candidate configuration is as follows: 


FCS(X) = 1— min {MF,(X)} (7) 
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Furthermore, each candidate configuration must be satisfied the following constraints: The radial 
topology of the EDS must be maintained and all of nodes must be served and all of nodes’ voltage must be in 
permissible ranges as well as no branch is over load. 


3. ICSA FOR THE MULTI OBJECTIVE NR PROBLEM 
The application of ICSA for the multi-objective NR is described in more details as follows: 
Step 1: Initialization 

At the beginning, the population of network configuration is generated as follows: 


X,= fr[X10,j +1. (Xnij = Xioi)| _ 


Where X; is the candidate configuration ith with i=1, 2, ..., N. N is population size. fr is a round function. r; 
is arandom number in [0, 1]. Xnjj; and Xo; are the upper and lower bounds of the switch jth with j=1, 2, ..., 
D. D is the number of opened switches in the EDS. 

From the created population, the fitness function shown in (7) is evaluated and the current best 
solution Xgbest is determined. 
Step 2: Creation of new population based on Lévy flight 

In this step, a new population network configuration is created by the Lévy flight as follows: 


XpPew = fr[Xbest; +1, x Axpe”] ©) 


Where, R2 is a random number in [0, 1]. AX}"°”is an increased value determined by the Lévy movement [20]. 
X best; is the so far best solution of the configuration ith in the population. 

Similar to step 1, the fitness function is evaluated. Then the current population is updated if the new 
individual is better than the corresponding once in the current population and the current best solution X goes: 
is updated. 

Step 3: Creation of new population based on alien eggs discovery 
In this step, there is also a new population generated as follows: 


Xj’ = fr[Xbest; + K(i,: )@AX/?™] ve 


Where AX}"*” is an increased value determined by the random movement [20]. K is an matrix with the value 
of elements of 0 or | controlled by the P, > 0 factor: 


K =rand (N,D) > P, (11) 


From the new population, the fitness function is evaluated. Then the current population is updated if 
the new individual is better than the corresponding once in the current population and the current best 
solution X pes; is updated. 

Step 4: Creation of new solutions around the current best solution 

In this step, there are D new solutions is created by modified one by one variable of the current best 
solution X goes. In which, the value of the modified control variable is controlled by the distance d from 
the current to new variable [27]. In ICSA, the value of d is set to the dina and dinin, So there are 2.D new 
solutions are created around the Xgbest. 


Xnew, = ft[Ky5.4.13 + Xgnese(O)] with j = 1,2,...,D (12) 


Where Kj.) is a vector with all elements of 0 except the element jth set to 1. R3 is arandom number in [0, 1]. 
In addition, to exploit more the search space around the X¢gpes, the Ky-) is also replaced by 
the random vector Krj.) which has more than one element set to 1. Therefore, there are also 2.D new 
solutions are created around the Xgpes: by modified more than one control variable of the Xgbest. From the 4.D 
new solutions, the fitness function is evaluated and the current best solution X ges: is updated. 
Step 5: Stopping the ICSA for the NR problem 
The searching processes of new configuration will be stopped as the number of evaluation 
the fitness function is reached to the pre-set value (MFE). In addition, because of taking long time for 
calculating the fitness function, the number of iterations called (UteFnot-improvea) that the fitness function value 
is not be improved is also considered as another stopping condition for ICSA. If two above conditions are not 
satisfied the ICSA will move to step 2 for searching elsewhere ICSA will be stopped and the best 
configuration Xgpes; is reported. 
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4. RESULTS AND DISCUSSION 

In this section, there are two illustrative system is used to reconfigure. The first system is used to 
validate the effectiveness of ICSA while the second system which is a practical system is used to show 
the application potential of ICSA for the practical systems. 


4.1. The 70-node system 

To evaluate effectiveness of ICSA for the multi-objective NR problem, the 70-node system 
consisting of 4 feeders, 78 branches and 11 opened switches shown in Figure 1 is used to reconfigure. 
The data of the system are referenced from [31]. In the initial configuration with opened switches of {69, 70, 
71, 72, 73, 74, 75, 76, 77, 78, 79}, the power loss of the system is 227.53 kW. The proposed ICSA method 
is implemented in Matlab software on the personal computer Intel(R) Core(TM) i5-2430 CPU @2.4GHz with 
4GB RAM. In addition, the runner root algorithm (RRA) is also run on the same computer to compare with 
ICSA. Due to successful application of RRA for the multi objective NR problem [32], RRA is selected for 
comparison. For ICSA, the Pz in (11), the dingy and dinin in (12) are set to 0.2, 4 and 2 respectively. For RRA, 
all parameters are selected similar to [32]. Some parameters selected the same for both ICSA and RRA 
are that the population size N=30, the stopping condition MFE=30000 and Iternot-improved=200. 





Figure 1. 70-node test system 


The performance of the ICSA is presented in Table 1. By using ICSA, the opened switches obtained 
in the optimal configuration are {46, 51, 67, 70, 71, 73, 75, 76, 77, 78, and 79}. After reconfiguration, 
the power loss has been reduced from 227.5625 to 205.6133 kW, the minimum voltage has also been 
improved from 0.9052 to 0.9283 p.u. Furthermore, LBI is reduced from 0.0083 to 0.0069 and the maximum 
load carrying coefficient is decreased from 0.3471 to 0.2890. Specially, to achieve the aforementioned 
results, it only takes 6 changing switches. The voltage and current profiles illustrated in Figure 2 shows that 
most of node voltages have been enhanced and the load balance among branches is more improved after 
performing reconfiguration by ICSA. 

Regarding the effectiveness of ICSA for multi-objective NR problems, Table 1 shows that although 
both ICSA and RRA have determined the optimal configuration, it is clear that the effectiveness of ICSA 
is superior to RRA. Specifically, in 50 independent runs, ICSA has found the optimal configuration in 24 
runs, while the RRA has reached the optimal solution only in 3 runs. These result leads to poor quality of 
obtained solution from RRA compared to ICSA. The maximum, mean and standard deviations (STD) of 
the adaptive function obtained from ICSA are much smaller than those of RRA. They are 0.5261, 0.2857 
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and 0.0427 respectively for ICSA while their values for RRA are 0.6951, 0.3870 and 0.1089. In addition, 
the average number of converged iterations of ICSA is also less than that of RRA. The maximum, minimum 
and mean convergence characteristics of both methods in 50 runs are shown in Figure 3. The figure 
demonstrates that the corresponding convergence curves of ICSA converge to smaller values than RRA. In 
addition, the mean convergence curve of ICSA is much closed to the minimum convergence curve. 
This shows the stability of ICSA for the multi-objective NR problem. 


Table 1. The comparisons among ICSA with RRA for the 70-node system 



























































Item Initial ICSA RRA 
Optimal configuration 69, 70, 71, 72, 73, 74, 46, 51, 67, 70, 71, 73, 46, 51, 67, 70, 71, 73, 
75, 76, 77, 78, 79 75, 76, 77, 78, 79 75, 76, 77, 78, 79 

Ploss (kW) 227.5625 205.6133 205.6133 

Vmmin (pu) (node) 0.9052 0.9283 (29) 0.9283 (29) 

LBI 0.0083 0.0069 0.0069 

Max(I/Irate) 0.3471 0.2890 0.2890 

LBF 0.1107 0.0086 0.0086 

‘ [1.0991, 1.1992, 1.7853, [1.2863, 1.4449, 1.4735, —_[1.2863, 1.4449, 1.4735, 

ean) 1.6375] 1.4865] 1.4865] 

NSW 0 6 6 

Number of runs found out the optimal _ 24/50 3/50 

solution/total runs 

Max of fitness - 0.5261 0.6951 

Min of fitness 1 0.2727 0.2727 

Mean of fitness - 0.2857 0.3870 

STD of fitness - 0.0427 0.1089 

The number of convergence iterations - 60 118 

CPU times (second) - 162.0978 97.7875 
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Figure 2. Voltage and current profile of the 70-node system before and after reconfiguration; (a) voltage 
profile, (b) current profile 
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Figure 3. Convergence characters of ICSA and RRA for the 70-node system 
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The comparisons among ICSA with other method in literature are shown in Table 2 also lead to 
the better obtained solution. Although the indexes of Vinin, LBI, LBF and NSW obtained from Multi-RRA [32] 
is nearly the same with ICSA, power loss obtained from Multi-RRA is still higher than that of ICSA. 
Compared to hybrid PSO (HPSO) [33], combination of heuristic rules and fuzzy multi objective approach 
(FMA) [31] and combination of honey bee mating optimization-discrete PSO (DPSO-HBMO) [34], all of 
indexes are similar to those of ICSA however it only takes 6 switching numbers compared to 12 of compared 
methods. In addition the balance of among memberships functions of ICSA compared to other methods 
shown in Table 3 demonstrated that the balance of among membership functions of ICSA is better than to 
others. Variance (var) of the membership function values obtained by ICSA is the smallest among 
the compared methods. 


Table 2. The comparisons among ICSA with other methods 

















Optimal opened AP Vinin Max 
Method So tchas (kW) (p.u.) LBI (Wtrate) LBF LBF; NSW 
46, 51, 67, 70, 71, 73, [1.2863, 1.4449, 
ICSA 75, 76, 77, 78, 79 205.61 0.9283 0.0069 0.2890 0.0086 1.4735, 1.4865] 6 
: 69, 37, 71, 67, 73, 45, [1.2807, 1.4449, 
Multi-RRA [32] 75,16 77,78, 79 208.31 0.9283 0.0064 0.2929 0.0098 1.4625, 1.5069] 6 
14, 79, 71, 39, 51, 28, [1.2871], 1.4471, 
HPSO [33] 73, 67, 46, 76, 70 205.07 0.9273 0.0066 0.2890 0.0085 1.4733, 1.4865] 12 
14, 79, 71, 39, 51, 28, [1.2871], 1.4471, 
FMA [31] 73, 67, 46, 76, 70 205.07 0.9273 0.0066 0.2890 0.0085 1.4733, 1.4865] 12 
DPSO-HBMO 14, 79, 71, 39, 51, 28, [1.2871], 1.4471, 
[34] 73, 67, 46, 76, 70 205.07 0.9273 0.0066 0.2890 0.0085 1.4733, 1.4865] 12 
Table 3. Balance among memberships functions of ICSA compared to other methods 
Method MFoPioss MFa-vmin MF 51 MEnsw MF pF var (MFi) 
ICSA 0.83911 0.83863 0.73684 0.72727 0.92455 0.0067 
Multi-RRA [32] 0.7359 0.8384 0.7182 0.72727 0.9139 0.0074 
HPSO [33] 0.8601 0.8016 0.6499 0.4545 0.9252 0.0355 
FMA [31] 0.8601 0.8016 0.6499 0.4545 0.9252 0.0355 
DPSO-HBMO [34] 0.8601 0.8016 0.6499 0.4545 0.9252 0.0355 





4.2. The 83-node practical system 

The second system which is a practical system of the Taiwan Power Company. There are 83 nodes, 
11 feeders and 13 open switches. The single line diagram and data of the system is referenced from [35]. At 
the initial radial configuration, power loss is 531.9924 kW and minimum voltage is 0.9285 p.u. Because of 
lack of rate current, the LBI membership is ignored. For ICSA, the population size and the stopping condition 
is chosen to 40 and 40000, respectively. Other control parameters of ICSA are selected similar to those of 
the 70-node system. 

The performance of the ICSA is presented in Table 4. By using ICSA, the opened switches obtained 
in the optimal configuration are {7, 34, 63, 84, 86, 87, 88, 89, 90, 91, 92, 93, and 95}. After reconfiguration, 
the power loss has been reduced from 531.9924 to 477.0592 kW, the minimum voltage has also been 
improved from 0.9285 to 0.9479 p.u. Furthermore, LBF is reduced from 1.4418 to 0.7124. The maximum 
power difference among feeders has decreased from 3.2840 to 2.7044 MVA. Specially, to achieve 
the aforementioned results, it only takes 6 changing switches. The voltage profiles illustrated in Figure 4 
show that most of node voltages have been enhanced after performing reconfiguration by ICSA. 


Table 4. The application of ICSA for the practical 83-node system 








Item Initial ICSA 
Open switches 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96 7, 34, 63, 84, 86, 87, 88, 89, 90, 91, 92, 93, 95 
Ploss (kW) 531.9924 477.0592 
Vmnin (pu) (node) 0.9285 (10) 0.9479 (84) 
LBF 1.4418 0.7124 
LBFi (MW) [4.4317, 3.3764, 4.5232, 2.8043, 4.6394, 1.3554, [2.5865, 3.3764, 4.5232, 2.8043, 2.9195, 3.0325 
3.3428, 1.8486, 3.5155, 1.9114, 4.6158] 4.1470, 2.8228, 3.5155, 1.9114, 4.6158] 
NSW 0 6 
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Figure 4. Voltage profile of the 83-node practical system 


5. CONCLUSION 


In this paper, the application of ICSA for the multi-objective NR problem has been successfully 
implemented for minimizing power loss, load balancing among the branches and among the feeders, node 
voltage deviation and switching operation numbers. The max-min method is used to selecte best compromise 
solution. The proposed method has been evaluated on the complex 70-node system. The numerical results 
validate that ICSA outperforms to RRA in terms of success rate, the obtained optimal solution quality and the 
number of converged iterations and other methods in literature in term of the quality of membership function 
values. In addition, ICSA is also applied for the practical system and the obtained results is more effective 


than the initial configuration. Therefore, ICSA can be the efficient and promised for solving the NR 
problems. 
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